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Higgs boson production via weak boson fusion at the CERN Large 
Hadron Collider has the capability to determine the dominant CP na- 
ture of a Higgs boson, via the tensor structure of its coupling to weak 
bosons. This information is contained in the azimuthal angle distri- 
bution of the two outgoing forward tagging jets. The technique is 
in dependent of both the Higgs boson mass and the observed decay 
^faannel. 

o : 

The CERN Large Hadron Collider (LHC) is generally re- 
garded as a tool that can guarantee direct observation of a Higgs 
jSbson, the remnant of the mechanism believed responsible for 
^electroweak symmetry breaking and fermion mass generation, 
£pnd the last unobserved element of the Standard Model (SM) of 
^ejementary particle physics. Furthermore, the LHC promises 
complete coverage of Higgs decay scenarios [|lj^], includ- 
general MSSM parameterizations [[l|,^||, and even invisible 
, Biggs decays [Q]. This capability has been greatly enhanced 
^jcently by the addition of the weak boson fusion (WBF) pro- 
(vdjiction channel to the search strategies |Q^J^]. While being 
l/eXtremely useful at the LHC, WBF has too low a rate and is too 
QJJmilar to background processes at the Fermilab Tevatron [R]. 
^-^ 0bservation of a resonance in some expected decay channel 
~-4s, however, only the beginning of Higgs physics. Continuing 
TSfforts will include the search for more than one Higgs boson, 
jtrpredicted e.g. by two-Higgs doublet models, of which the 
{M$SM [||,||] is a subset. At least as important is the detailed 
%Judy of the properties of the Higgs-like resonance, not only at 
'TSlfuture Linear Collider jTo| ] but also at the LHC: determination 
all the quantum numbers and couplings of the state. These 
kiflclude the gauge, Yukawa and self-couplings as well as the 
^Jiarge, color, spin, and CP quantum numbers. While charge 
Cahd color identification is straightforward and a technique has 
been proposed for the gauge and Yukawa coupling determina- 
tions [fyl], the LHC has considerable difficulty in practice to 
determine the Higgs CP transformation properties for interme- 
diate Higgs masses [ |T2[ ] via a weak boson coupling, and no 
technique has yet been proposed to identify the tensor structure 
of the Higgs-weak boson vertex in the intermediate mass range. 
The methods ofRef. [[DJ may be useful, but only for very light 
Higgs masses. Furthermore, this method does not involve the 
weak boson vertices at all. 

In this letter we propose a technique which achieves the CP 
measurement goal via a study of WBF events. WBF Higgs 
production, while not the largest cross section at the LHC, is 
useful because of its characteristic kinematical structure, in- 
volving two forward tagging jets and central Higgs decay prod- 
ucts, which allows one to isolate the signal in a low background 



environment. The angular distribution of the two tagging jets 
carries unambiguous information on the CP properties of the 
Higgs couplings to weak bosons which is independent of the 
Higgs decay channel observed. 

As a theoretical framework we consider two possible ways 
to couple a spin zero field to two gauge bosons via higher di- 
mensional operators. In a gauge invariant dimension six (D6) 
Lagrangian, the terms 
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lead to anomalous couplings between the Higgs-type scalar and 
two charged gauge bosons [JfJI]. The scales A c and A Q set the 
coupling strength of CP even and CP odd scalars, respectively. 
The Feynman rules can be read off the dimension five (D5) 
operators that result when $ is given a physical field expansion: 
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and similarly for the Z boson. The two scales are related via 
I/A5 = g 2 v/K\. Since we assume SU(2) invariance and do 
not consider additional D6 operators like B^B^ , the 
WWH and ZZH couplings are related by the same cos 2 6w 
factor as in the SM. 

In principle one would have to introduce a form factor to en- 
sure the unitarity of scattering amplitudes involving these op- 
erators. However, we have checked that at the LHC the typical 
Pt of the tagging jets, for WBF processes generated by the D5 
operators, remains comparable to the SM case and is well be- 
low the scale A, which we assume to be of order a few hundred 
GeV or above. Thus, form-factor effects would remain small 
in a more complete treatment and they would not distort the 
angular distributions to be discussed below. 

The analog of the CP even D5 operators is present in Higgs 
production through gluon fusion, as HG^G^" ', and gives 
an excellent approximation for the ggH coupling induced by 
heavy quark (and squark) loops. In the low energy limit the 
D5 operators also appear in the one-loop WWH coupling, 
but their size is suppressed by a factor aw /tt ~ 10 -2 and 
hence not observable at the LHC, as we will see later. Another 
source would be a Higgs-like top-pion that is a general feature 
of topcolor models [[l5|] and which couples to weak bosons like 
HWf v W~ with a coefficient that is considerably larger than 
in the SM and is expected to lead to observable rates of produc- 
tion in the WBF channel. 
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For a true Higgs boson the WWH and ZZH couplings orig- 
inate from the kinetic energy term of the symmetry breaking 
field, (_D A1 $)t(Z)' i $), which mediates couplings proportional 
to the metric tensor. This tensor structure is not gauge-invariant 
by itself and identifies the Higgs field as the remnant of spon- 
taneous symmetry breaking. It is thus crucial to distinguish it 
from the effective couplings derived from Eq. ([!]). Since the 
partons in the WBF processes 

VP ~ * q<l'H — > qq'TT, qq'WW, qq'^fj (3) 

are approximately massless, the production cross section is 
proportional to the Higgs-weak boson coupling squared. Re- 
placing the g^" coupling with a higher dimensional coupling 
changes the kinematical structure of the final state scattered 
quarks. 

To illustrate this we consider leptonic final states in H — > rr 
decays as in Ref . [^| . We emphasize the H — > rr decay chan- 
nel because it is resilient to modifications of the Higgs sector as 
encountered in the MSSM: a luminosity of 40 fb -1 guarantees 
coverage of the entire (m^-tan (3) plane after combining the 
leptonic and semileptonic decay channels of the tau pair [||]. 
The basic set of cuts on the outgoing partons consists of 

PTj > 20 GeV AR j:j > 0.6 |^| < 4.5 

in addition to the separation and acceptance cuts for the de- 
cay leptons, which we don't discuss here. (Further cuts on the 
invariant mass of the tagging jets and the tau pair decay kine- 
matics are necessary to extract the signal. These details and the 
final step of reconstructing the tau pair invariant mass are cur- 
rently under study by various CMS and ATLAS groups, with 
very encouraging results flip]].) In the parton level analysis we 
are left with a cross section of a ~ 0.5 fb for a 120 GeV SM 
Higgs boson, leading to S/B = 2.7/1 and a Gaussian signif- 
icance (TGauss = 6.8 for 60 fb -1 of data [Q]. The two largest 
backgrounds are QCD and electroweak rrjj production, which 
together are < 30% of the signal cross section after cuts. The 
other backgrounds, including H — > WW and tt+ jets, are of 
minor importance and can safely be neglected in the following 
qualitative analysis. 

1. Let us first assume that a Higgs-like scalar signal is found 
at the LHC in this channel at the expected SM rate. We must 
experimentally distinguish a SM g^-type coupling from the 
tensor structures implied by the D5 operators of Eq. (||). A 
SM rate induced by one of the D5 operators requires a scale 
A 5 w 480 GeV (A 6 w 220 GeV). A particularly interesting 
kinematic variable is the azimuthal angle A^jj between the two 
tagging jets. For forward scattering, which is dominant due to 
the M^-propagator factors, the remaining SM matrix element 
squared for qq — > qqH is proportional to s rrfc -, where nijj is 
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Figure 1. Azimuthal angle distribution between the two tagging 
jets for the signal and dominant rr backgrounds, itih = 120 GeV. 
Cross sections for the D5 operators correspond to A5 = 480 GeV, 
which reproduces the SM cross section, after cuts as in Eq.(Q) and 
Ref. The expected SM background is added to all three Higgs 
curves. 



the invariant mass of the two tagging jets. This leads to an es- 
sentially flat azimuthal angle distribution between the two jets, 
as shown in Figs, [j] and |f In the H — > rr case, a slight bias 
toward small angles is introduced by selection cuts, which re- 
quire a substantial transverse momentum for the Higgs boson. 
The major backgrounds, Zjj production with Z — > rr, possess 
mostly back-to-back tagging jets. 

For the CP odd D5 operator, the shape of the distribution 
follows from the presence of the Levi-Civita tensor in the cou- 
pling: it gives a nonzero result only if there are four indepen- 
dent momenta in the process (here, the four external parton mo- 
menta). For planar events, i.e. for tagging jets which are back- 
to-back or collinear in the transverse plane, the matrix element 
vanishes. 

The CP even operator given in Eq.(||) develops a special fea- 
ture for forward tagging jets. In the limit of |pi tag ' ) | 3> \pxaP \ 
and small energy loss of the two scattered quarks, we can ap- 
proximate the matrix element by 

Mc,5 oc — !— Jf J2 [g^{qi ■ 92) - qxuqin] 

A e,5 

- t~ [4 4 - JfJi] p^ tagl) • p? ag2) (5) 

A c , 5 

where q t , Ji are the momenta and currents of the intermediate 
weak gauge bosons. For A0jj = tt/2 the last term vanishes, 
leading to an approximate zero in the distribution. From the 
three curves shown in Fig. [l] we conclude that the azimuthal 
angle distribution is a gold plated observable for determining 
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Figure 2. Normalized distributions of the azimuthal angle be- 
tween the two tagging jets, for the H — > WW — > e/i^T signal at 
m,H = 160 GeV. Curves are for the SM and for single D5 operators 
as given in Eq. (^), after cuts as in Eq.(Q) and Ref. [Q]. 



the dominant CP nature and the tensor structure of the Higgs 
coupling. With 100 fb" 1 of data per experiment, the SM case 
can be distinguished from the CP even (CP odd) D5 couplings 
with a statistical power of ~ 5 (4.5) sigma, from the H — > rr 
channels [^[] alone. This observable is furthermore independent 
of the particular decay channel and Higgs mass range. We 
have explicitly checked the case of a 160 GeV Higgs boson 
decaying to W pairs and find exactly the same features, shown 
in Fig. [| Note, however, that in this case decay distributions 
will depend on the structure of the HWW vertex also. 



2. Let us now examine the following scenario: a Higgs 
candidate is found at the LHC with a predominantly Standard 
Model g^ v coupling. How sensitive will experiments be to any 
additional D5 contribution? 

For the CP odd D5 coupling we do not observe any inter- 
ference term between the Standard Model and the D5 matrix 
element. Although there is a non-zero contribution at the ma- 
trix element level, any hadron collider observable is averaged 
over charge conjugate processes since we cannot distinguish 
quark from antiquark jets. As a result, interference effects can- 
cel in typical hadronic differential cross sections. Using the 
azimuthal angle distribution will only marginally enhance the 
sensitivity to a small contribution of the CP odd Higgs coupling 
beyond what a measurement of the Higgs production cross sec- 
tion could give. 

In the case of a contribution from a CP-even D5 opera- 
tor, interference effects are important for the distortion of the 
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Figure 3. Azimuthal jet angle distribution for the SM and inter- 
ference with a CP even D5 coupling. The two curves for each sign of 
the operator correspond to values ij/<jsm = 0.04, 1.0. Error bars for 
the signal and the dominant backgrounds correspond to an integrated 
luminosity of 100 fb -1 per experiment, distributed over 6 bins, and 
are statistical only. 



distribution. 



All additional terms in the squared ampli- 
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tude \M\ = \M.sm + Mc.sY nave an approximate zero at 
A0jj = 7r/2, according to Eq.(||). Moreover, the dominant 
piece of the anomalous amplitude changes sign at this approx- 
imate zero which results in a sign change of the interference 
term at n/2. Fig. |3] shows that, dependent on the sign of the 
D5 operator, the maximum of the distribution is shifted to large 
or small angles A^jj . Results are shown for two different val- 
ues of the scale A 5 which are chosen such that the D5 operator 
alone, without a SM contribution, would produce a Higgs pro- 
duction cross section, a, which equals 0.04 (1.0) of the SM 
cross section, ctsm- While changes in cross sections of a few 
percent are most likely beyond the reach of any LHC count- 
ing experiment, we see that in the differential cross section the 
effect of D5 operators is quite significant [[17]]. 

To quantify this effect and at the same time minimize sys- 
tematic errors we define the asymmetry 
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<r(A^jj < tt/2) - <r(A^jj > tt/2) 
^(A^jj < tt/2) + ^(A^j > tt/2) 



(6) 



One major source of systematic uncertainty will be the gluon 
fusion induced H + 2 jets background, which in the large 
top mass limit is proportional to the CP even D5 operator 
HG ^ V G^ V . At the amplitude level, this operator induces the 
same azimuthal angle dependence of the two jets as the CP even 
operator of Eq. (||). However, since it contributes to H + 2 jets 
via ^-channel gluon (color octet) exchange, it cannot interfere 
with WBF This gluon fusion contribution can exceed O(10%) 
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Figure 4. Comparison of the sensitivity of a total cross section 
(counting) experiment and the azimuthal angle asymmetry, Eq.(^), to 
the presence of a CP even D5 coupling. The horizontal lines rep- 
resent one sigma statistical deviations from the SM value. The sec- 
ondary axes show the corresponding values of A5 and A6, as defined 
mEqs.(R§). 



of the signal after cuts [ |T8|] and is expected to have large higher 
order QCD corrections Jl9|]. The measurement of the absolute 
rate of WBF events would therefore be systematics limited, due 
to the unknown if-factor for the gluon fusion contamination. 
Assuming that this i^-factor does not vary with A<frjj, a full 
shape analysis of the azimuthal angle distribution allows to dis- 
tinguish this noninterfering gluon fusion background from an 
interfering D5 HWW coupling: the asymmetry is dominated 
by the interference terms. As mentioned before there is a loop 
induced WWH coupling, but it is expected to contribute with 
size a ~ (a/ir) 2 <tsm> beyond the reach of even a linear col- 
lider precision experiment [ |Tot l . 

In Fig. [| we compare the sensitivity to D5 couplings ex- 
pected from the total cross section and the azimuthal asym- 
metry, respectively. In the integrated cross section, interfer- 
ence effects between the SM coupling and the CP even 
D5 coupling largely cancel. With 100 fb -1 per experiment, 
a total cross section measurement at the LHC is sensitive (at 



the l-cr level, and considering statistical errors only) to A c .e < 
510 GeV (1/A > 0) or 290 GeV (1/A < 0). In contrast,' 
is a much more sensitive observable, and equally sensitive to 
positive and negative A. For both signs of the D6 coupling the 
reach in the leptonic rr channel is ~ 690 GeV, significantly 
better than the counting experiment. A rough estimate shows 
that, for a 120 GeV Higgs boson, the LHC will be sensitive 
to Ag ~ 1 TeV, after adding the statistics of both rr the 
WW [|], and the 77 [g] WBF channels. While this reach is 
not quite competitive with the linear collider analysis Jl0[], it is 
by no means certain that a linear collider will be built, and our 
work adds a dimension to the LHC Higgs analysis previously 
thought not possible. 

In summary, the weak boson fusion production process is 
not only a competitive discovery channel for an intermediate 
mass Higgs boson, it also offers the opportunity to unveil the 
structure of the Higgs field's coupling to gauge bosons. Using 
information obtained with generic weak boson fusion cuts for 
the intermediate-mass Higgs search, one can unambiguously 
determine the CP nature of a Higgs-like scalar: the azimuthal 
angle distribution between the tagging jets clearly distinguishes 
the Standard Model coupling from a typical loop induced 
CP even or CP odd coupling. In a search for dimension five op- 
erators which interfere with the SM HWW coupling, an asym- 
metry analysis of this azimuthal angle distribution improves the 
reach far beyond what is possible in a counting experiment, in- 
cluding the determination of the sign of the additional coupling. 



ACKNOWLEDGMENTS 

We want to thank T. Han and O. Eboli for inspiring dis- 
cussions. This research was supported in part by the Univer- 
sity of Wisconsin Research Committee with funds granted by 
the Wisconsin Alumni Research Foundation and in part by the 
U. S. Department of Energy under Contract No. DE-FG02- 
95ER40896. 



[1] Z. Kunszt and F. Zwirner, Nucl.Phys. B385 (1992) 3; AT- 
LAS TDR, report CERN/LHCC/99-15 (1999); CMS TP, re- 
port CERN/LHCC/94-38 (1994); M. Spira, Fortschr.Phys. 46 

(1998) 203 and references therein. 

[2] V. Barger, G. Bhattacharya, T. Han, and B.A. Kniehl, Phys.Rev. 
D43 (1991) 779; M. Dittmar and H. Dreiner, Phys.Rev. D55 
(1997) 167. 

[3] D. Rainwater, D. Zeppenfeld, and K. Hagiwara, Phys.Rev. D59 

(1999) 014037; T. Plehn, D. Rainwater, and D. Zeppenfeld, 
Phys.Lett. B454 (1999) 297; Phys.Rev. D61 (2000) 093005. 

[4] O.J.P Eboli and D. Zeppenfeld, Phys.Lett. B495 (2000) 147. 
[5] D. Rainwater and D. Zeppenfeld, JHEP 9712 (1997) 5. 



4 



[6] D. Rainwater and D. Zeppenfeld, Phys.Rev. D60 (1999) 
113004, erratum ibid. D61 (2000) 099901; N. Kauer, T. Plehn, 
D. Rainwater, an d D. Zeppenfeld Phys.Lett. B503 (2001) 113. 

[7] M. Carena et at., hep- ph/0010334 

[8] J.F. Gunion, Phys.Lett. B322 (1994) 125; V. Barger, 
RJ. Phillips, and D.P. Roy, Phys.Lett. B324 (1994) 236; 
K.A. Assamagan and Y. Coadou, ATL-PHYS-2000-031. 

[9] A. Djouadi, W. Kilian, M.M. Muhlleitner, and P.M. Zer- 
was, Eur.Phys.J. C10 (1999) 27; Eur.Phys.J. C10 (1999) 
45; T. Plehn, M. Spira, and P.M. Zerwas, Nucl.Phys. 
B479 (1996) 46; erratum ibid. B531 (1998) 655; S. Dawson, 
S. Dittmaier and M. Spira, Phys.Rev. D58 (1998) 115012. 
[10] TESLA Technical Design Report; DJ. Miller, S.Y. Choi, 
B. Eberle, M.M. Muhlleitner, and P.M. Zerwas, Phys.Lett. 
B505 (2001) 149; K. Hagiwara and M.L. Stong, Z.Phys. C62 

(1994) 99; M. Kramer, J. Kiihn, M.L. Stong, and P.M. Zerwas, 
Z.Phys. C64 (1994) 21; K. Hagiwara, S. Ishihara, J. Kamoshita, 
and B.A. Kniehl, Eur.Phys.J. C14 (2000) 457; T. Han and 
J. Jiang, Phys.Rev. D63 (2001) 096007. 

[11] D. Zeppenfeld, R. Kinnunen, A. Nikitenko, and E. Richter-Was, 

Phys.Rev. D62 (2000) 013009. 
[12] The original proposals for CP measurement techniques may be 

found in: J.R. Dell'Aquila and C.A. Nelson, Phys.Rev. D33 

(1986) 80; Phys.Rev. D33 (1986) 93; C.A. Nelson, Phys.Rev. 

D37 (1988) 1220; and references therein. 
[13] J. F. Gunion and X. G. He, Phys. Rev. Lett. 76, 4468 (1996). 
[14] See e.g. W. Buchmuller and D. Wyler, Nucl.Phys. B268 

(1986) 621; K. Hagiwara, R. Szalapski and D. Zeppenfeld, 

Phys.Lett. B318 (19 93) 155. 

[15] See e.g. G. Burdman, frep-ph/961 1265] ; K. Lane and E. Eichten, 

Phys.Lett. B352 (1995) 382; C.T Hill, Phys.Lett. B345 

(1995) 483; and references therein. 

[16] See e.g. the talks by K. Jakobs and A. Nikitenko at the "Work- 
shop on the Future of Higgs Physics", May 3-5, 2001, Fermilab. 

[17] M.C. Gonzales-Garcia, Int.J.Mod.Phys. A14 (1999) 3121; 
O.J.P. Eboli, M.C. Gonzalez-Garcia, S.M.Lietti, and S.F. Novaes, 
Phys.Lett. B478 (2000) 199. 

[18] V. Del Duca, W.B. Kilgo re, C. Oleari, C. Schmidt, and D. Zep- 
penfeld, hep-ph/0105129| . 

[19] M. Spira, A. Djouadi, D. Graudenz and P.M. Zerwas, 
Nucl.Phys. B453 (1995) 17; M. Kramer, E. Laenen and 
M. Spira, Nucl.Phys. B511 (1998) 523; S. Catani, D. de Flo- 
rian and M. Graz zini, JHEP 0105 (2001) 25; R.V. Harlander 



and W.B. Kilgore, iep-ph/0 102241 



5 



